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We theoretically calculate the average fraction of frozen particles in rectangular systems of arbi- 
trary dimensions for the Kob- Andersen and Fredrickson- Andersen kinetically-constrained models. 
We find the aspect ratio of the rectangle's length to width, which distinguishes short, square-like 
rectangles from long, tunnel-like rectangles, and show how changing it can effect the jamming tran- 
sition. We find how the critical vacancy density converges to zero in infinite systems for different 
aspect ratios: for long and wide channels it decreases algebraically v c ~ W -1 ^ 2 with the system's 
width W, while in square systems it decreases logarithmically v c ~ 1/lnL with length L. Although 
derived for asymptotically wide rectangles, our analytical results agree with numerical data for 
systems as small as W ~ 10. 

PACS numbers: 45.70.-n,64.60.an,64.70.Q- 



I. INTRODUCTION 

Increasing the density of particles in granular matter 
causes them to undergo a transition from an unjammcd 
state, where the particles can move relatively freely, to 
a jammed state, where almost none of the particles can 
move [l| . Systems of interest in nature and in industrial 
applications typically have complicated geometries which 
strongly affect jamming in them @-@|, and it is thus im- 
portant to understand how does confinement influence 
the jamming of granular matter. Here we investigate the 
effects of confinement on the jamming transition, and in 
particular test how does the shape of containers deter- 
mine how they jam. Most theoretical work so far was 
done on square systems [r34ITl|. 

There are numerous laboratory experiments that deal 
with non-square systems [134151 . For example, Daniels 
and Behringer conducted an experiment on polypropy- 
lene spheres in an annulus [l6j |. which is large enough to 
be considered a rectangle with infinite length and finite 
width. A different experiment by Bi et al. j!7| consists of 
shearing a square system such that it becomes a rectangle 
with the same area and particle density as the original 
square. In this paper we study the effects of confinement 
on jamming by studying these phenomena in kinetically- 
constrained models in rectangular domains. 

The essence of jamming is captured by the various 
kinetically-constrained models [l8l43ll . F or such simple 
models, and for other related models I32443I I. it is useful 
to have exact solutions. 



A. Kinetically-Constrained Models 

Two types of kinetically-constrained models which de- 
scribe granular and glassy materials are spin-facilitated 
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models and lattice-gas models [18|, ll9j ■ In both types of 
models, the system is represented by a grid, such that 
each site on the grid can have one of two values, 1 or 0. 

In the lattice-gas models, a site with a value of 1 rep- 
resents a single particle and a site with a value of rep- 
resents a vacant region. In each time step, one particle 
and one of the possible directions are chosen randomly 
with equal probabilities. The chosen particle attempts 
to move in that direction, and if the kinetic constraint 
allows the move, the particle moves to the neighboring 
site in the chosen direction. For a given initial configura- 
tion, some of the particles can move from the start, and 
some can move only after (many) other particles have 
moved and cleared the way for them. There may also be 
particles that will never move, no matter how the other 
particles in the system move. Those that will never move 
are called permanently frozen, and those that can move 
eventually are called unfrozen. See Fig. [T]for an example. 

In the case of spin-facilitated models, a site with a value 
of represents a region of low density and high mobility, 
and a site with a value of 1 represents a region of high 
density and low mobility. Note that this notation is dif- 
ferent than the common notation (0 for a region of high 
density and 1 for a region of low density), but we use 
this definition in order to deal with both lattice-gas and 
spin-facilitated models simultaneously. In each time step 
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FIG. 1: The difference between frozen and unfrozen par- 
ticles in the Kob- Andersen model: White square are vacan- 
cies. Light-gray particles can move in this initial configura- 
tion. Dark-gray particles cannot move now, but after some 
other particle(s) move, they too are mobile. Black particles 
are permanently frozen and will never move. 



of the dynamics, one of the sites is chosen randomly, and 
changes its value at a temperature-dependent probability 
if it has enough neighbors with low density (i.e., a value 
of 0), with the exact geometric criteria depending on the 
specific model. In the limit of zero temperature, the only 
allowed changes are from high density to low density, i.e. 
from a value of 1 to 0. For a given initial configuration, 
there is a possibility that even after an infinite number of 
time steps, some sites will still have a value of 1. These 
permanently frozen sites represent the backbone of the 
system which will never change. 

In order to find the fraction of permanently frozen 
particles, one can use the bootstrap method, which it- 
eratively removes mobile particles, until none of the re- 
maining particles can move. Again, we have a backbone 
of sites which will never change. This algorithm is obvi- 
ously valid for spin-facilitated models, but also for lattice 
gas models, since the criteria for the mobility of particles 
is local, and removing a mobile particle is effectively the 
same as moving it far enough from its neighbors. Since 
the algorithm for finding the backbone of both types of 
models is similar (but not identical), we will use the same 
language to describe both models, and choose the lan- 
guage of lattice-gas models. This means, for example, 
that whenever we speak of "vacancies" it should be in- 
terpreted as "sites with value 0" or "low-density regions" 
in the context of spin-facilitated models. For brevity, we 
will also use the term frozen particles interchangeably 
with permanently frozen particles. 

We consider a two-dimensional rectangle, represented 
by a square lattice, such that each site either contains 
one particle or is vacant. The rectangle has L sites in 
the horizontal direction, and W sites in the vertical di- 
rection, such that L > W . In our numerical simulations 
we used hard-wall and periodic boundary conditions in 
both directions, but most of our analytical approxima- 
tion ignores the boundary conditions. For rectangles of 
infinite length, hard-wall boundary conditions simulate 
particles inside a two dimensional channel, and periodic 
boundaries simulate particles on the surface of a cylin- 
drical tube. 

For the lattice-gas model we use the Kob-Andersen 
(KA) model [Hj], such that a particle can move if it 
has at least two neighboring vacancies before and af- 
ter the move. For the spin-facilitated model we use the 
Fredrickson- Andersen (FA) model (45|, such that a site 
can change its state if it has at least two neighboring sites 
with a value of 0. We could have chosen a different num- 
ber of neighbors needed for movement, but on the square 
lattice the only interesting effects occur at two neigh- 
bors. If only one neighbor is needed for movement then 
all the particles are movable as long as there is at least 
one vacancy in the FA model or two adjacent vacancies 
in the KA model. If three neighbors are needed for move- 
ment then any closed loop is frozen, even a 2 x 2 block, 
which means that almost all the particles in the system 
are frozen. Mathematically speaking, the KA model and 
the FA model are very similar to each other. Also, a 



mobile particle in the KA model is necessarily mobile in 
the FA model, and a frozen particle in the FA model is 
necessarily frozen in the KA model, thus the fraction of 
frozen particles in the KA model is larger than (or at 
least equal to) the fraction of frozen particles in the FA 
model. 

Toninelli, Biroli and Fisher showed [22| that for an 
infinite system in the KA model, none of the particles 
are permanently frozen as long as the lattice is not com- 
pletely full with particles, which automatically means 
that this is also the case in the FA model. In this paper 
we study how many particles, on average, are perma- 
nently frozen for a given particle density, p, and given 
rectangle dimensions W x L. 



B. Finite Size Effects 

Numerical simulations done on square systems 
(L = W, in our notations) showed that the fraction of 
permanently frozen particles, npF , rises rapidly from 
to 1 at a certain critical density, p c , which increases with 
system size. Holroyd [13] theoretically analyzed jamming 
in this context using the notion of critical droplets, which 
are small unjammed regions which facilitate movement 
throughout the system. He showed that for very large 
squares in the FA model the relation between the critical 
density and the system size is 

square -i _ 

Pc lni' 



Toninelli et. al. showed [22j that this value of A is also 
true for squares in the KA model. However, this result 
is only true for asymptotically large L. 
We can define an effective A 

A— e5 (L) = [l-p c (L)]lnL, (2) 

which converges to A as the system size increases. For 
systems of size L « 10 2 ~ 10 5 it was found [48| that 
^squares ^ q ^ f or a r] simulated sizes. This contradic- 
tion was resolved by Holroyd's proof [49| that the conver- 
gence of \ s e q j U j°' res to A is very slow and may be apparent 
only at systems of size L ps 10 20 , beyond the capabilities 
of modern computers, and beyond the range of physi- 
cal realization (L = 10 20 implies a system of L 2 = 10 40 
particles) . 

Holroyd's analysis considered only the size of the sys- 
tem, and not its shape. For long rectangular domains, 
this method may not be used. Instead, we find that rect- 
angular systems may be divided into independent sec- 
tions, in the sense that jamming in one section does not 
depend on the internal configuration within its neighbor- 
ing sections. Within each section, Holroyd's notion of 
critical droplets may be used. 
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C. Outline 

In this paper we show how not only the size of a system 
influences jamming in it, but also its shape. We demon- 
strate this, first by considering the limit of very large sys- 
tems. When considering a square system of size L x L, 
and taking the limit L — > oo, one finds that the critical 
vacancy density, v c = I — p c , scales as v c ~ 1/lnL. We 
find that when the system's width W is fixed and the 
length is taken to infinity L — oo, the critical density 
scales as v c ~ 1/y/W when W — > oo . 

The second scenario we consider is of a system of fixed 
area, for which we change the aspect ratio between the 
width and the length of the system. We find that stretch- 
ing the system causes it to jam, and relate this result to 
recent experiments of sheared granular matter. 

The paper is organized as follows. In Section [XT] of this 
paper we derive an approximate analytical expression for 
the fraction of frozen particles, npp, for a rectangular 
system of arbitrary dimensions W x L. In Section ITO1 we 
deal with large systems, R^X>1, and use our approx- 
imation to find the critical density at which the system 
goes from jammed to unjammed and the width of this 
transition. We find that the system can be considered 
infinite if its length is longer than the average section 
length, or equivalently In L ^> V4XW. In Section |IV] we 
deal with narrow systems (small W), and improve the 
approximation derived in Section [XT] We even derive an 
exact result for the case of very narrow systems (L = oo 
and W = 1,2). In Section [V] we investigate the internal 
structures in the system. The Appendices contain the 
derivation of the lengthy expressions used in the analyt- 
ical approximation. 

II. CRITICAL DROPLETS AND DIVISION 
INTO SECTIONS 

A. Critical Droplets 

Holroyd showed that in a large enough square there is 
a probability of approximately e~ 2X l v that a particle is 
part of a critical droplet, where A is given in Eq. ((lb). 
Hence, the total number of critical droplets in a rectangle 
of size L x W is WLe~ 2X / v , where v = 1 — p is the 
vacancy density. The expression for the critical density 
in square systems, Eq. (QJi), is derived by setting W = L. 
Such a critical droplet can cause the entire system to be 
unfrozen, thus the critical vacancy density is when the 
number of critical droplets is finite, since below (above) 
that critical density the number of critical droplets is very 
small (large) when L is taken to infinity. 

The fraction of frozen particles in the FA and KA mod- 
els is obviously different, due to the different kinetic con- 
straints, but that difference is small. The reason that the 
densities of frozen particles for both the FA and the KA 
models are almost the same can be seen from the con- 
struction of Holroyd's proof. Holroyd considered small 




FIG. 2: (color online) Division into sections, represented by 
different colors: Columns 1 — 2, 3 — 9, and 10 — 15. The strips 
are columns 1 — 2, 7 — 9, and 14 — 15. 

critical droplets, which are unfrozen, and checked how 
they can be expanded to "unfreeze" the entire system. 
The criterion for the expansion of these droplets is the 
same in both models, and the only difference is in what 
constitutes a small critical droplet. For example, the 

structure ^ ^ is unfrozen in the FA model but frozen 

in the KA model. For large enough systems, and evi- 
dently for small ones too, the effect of this difference is 
negligible. 

However, when L — > oo and W remains constant we 
cannot simply set 1 = W Le~ 2X / Vc to find the critical den- 
sity, since the solution to this equation is v c = for all W . 
Our resolution of this problem is obtained by dividing the 
long rectangle into finite sections, implementing the idea 
of critical droplets in each section, and finally averaging 
over all sections. Another approach, which yields the 
same results, is solving the equation 1 = W {I) e~ 2X / Vc , 
where (I) is the average section length (see below). 

B. Division into Sections 

A rectangular system may be divided into sections by 
noting that if there are two or more successive columns 
which are completely full, then all the particles in them 
are permanently frozen in both the KA and FA models 
and with either hard-wall or periodic boundary condi- 
tions. We call a pattern of m successive full columns a 
strip of size m, where m > 2. A single full column is not 
called a strip. These strips divide the rectangle into finite 
sections, such that the leftmost column on each section is 
the first not-full column after a strip, and the rightmost 
column is the final column of the next strip. Each section 
contains only one strip. For example, see Fig. [2J 

Using this division, we note that the particles within 
one section are independent from the particles in the 
other sections, i.e., the state of the particle (whether it 
is frozen or not) depends only on the structure within 
its section, and not on the configuration of neighboring 
sections. 

Each finite section of length l+m ending with a strip of 
size TO and with n + mW occupied sites, has many config- 
urations for the 77 particles in the I columns not occupied 
by the strip. We will denote these configurations with an 
index s. Since the probability of having a strip of length 
777 is independent of the probability for a certain config- 
uration in the rest of the strip, the probability of such a 
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configuration occurring is 

P(n,l,m,s) = p mW Q(n,l,s), 



<N> 



(3) 



where p mW is the probability of having a strip of length 
m containing mW particles, and Q(n,l, s) is the relative 
probability of configuration s in the region with / columns 
and n occupied sites between the strips, such that there 
are no two adjacent full columns. The reason we exclude 
the possibility of two adjacent full columns is to count 
each type of section only once, since two (or more) adja- 
cent full columns divide the section into smaller sections. 
The average fraction of frozen particles, npp, is the num- 
ber of frozen particles divided by the number of particles, 



npF 



(4) 



where N(n, I, m, s) — N(n, I, s) + mW is the total num- 
ber of permanently frozen particles in the section, with 
N(n,l,s) being the number of frozen particles in the I 
columns not occupied by the strip. The sum over I and 
m is such that I + m < L, and m is greater or equal to 
2, except in the following special cases: no strip in the 
entire rectangle (I = L,m = 0), and the entire section is 
full (/ = 0, m = L). Since the probability P appears both 
in the nominator and the denominator in Eq. ((3]) we need 
not worry about its normalization or the normalization 
of Q. However, we find that in the limit of infinite length 
the probability is normalized such that Y]P = 1. 



C. Our Approximation for Rectangular Systems 

In our case, we assume that the probability that a par- 
ticle in a section of length / + 1 is frozen is the probability 
that it is frozen in a section of length I times the prob- 
ability that the added column does not contain critical 
droplets, 



(N(l + 1)) = (N(l)) (l-e- 2X ^ 



w 



(5) 



where (..) is the average over all configurations. The 
solution to this recursion relation is 

(N(n,l)) = ne~ w , k = -W In (l - e" 2A/ ") . (6) 

This leads to very good agreement with results of numer- 
ical simulations, as shown in Fig. [3] 
Using this assumption in Eq. ((3]) yields 

_ E n ,i, m P mW Q(n,l) [mW + ne- kl ] 

where Q(n,l) — J2 s Q( n ' '> s )- ^ ne evaluation of these 
sums in closed form is given in Appendix A. We note 
here that the denominator and the first part of the nom- 
inator in Eq. ([7]) do not depend on our approximation 




1 



4 w 



FIG. 3: (color online) Average fraction of frozen particles 
in sections of length I in the KA model with hard- wall bound- 
ary conditions for different widths and densities: W = 4, 
p — 0.8 (blue squares), W = 10, p — 0.89 (purple circles), 
W = 20, p = 0.93 (yellow triangles), W = 40, p = 0.95 (green 
stars). Continuous lines are approximations with A taken 
from simulations of long rectangles (L = 200 W). 



relating the fraction of frozen particles with the number 
of critical droplets, and thus are not approximations but 
exact results. We further note that the ratio between the 
first part of the nominator and the denominator is the 
density of particles which are in the strips in both mod- 
els and with both boundary condition. Hence, we define 
the density of particles which are in the strips as 



''strip 



T, n ,i,mP mW Q(nJ) [mW + n] 



which in the limit of L — ¥ oo converges to 



n str ip(L -t oo) = p 



2W-1 



(2 



P 



in 



(8) 



(9) 



We find that even for W as small as 10 and for all L > W, 
the density of particles in the strips is very low, except in 
the region very near p = 1, where the fraction of frozen 
particles, npp, is almost unity. This means that for wide 
systems, the strips hardly contribute to the total fraction 
of frozen particles near the critical density and below it. 
The only role the strips play in this regime is dividing the 
system into sections, which are very long since there are 
few strips. Using the results in Eqs. (|A17I) and (|A20|1 . 
the fraction of frozen particles at L = oo can be written 
as 



npp(L 



where 



J2W—1 



iW-l 

p w ) + ^^N PF (p,W) 7 
(10) 



N PF (p, W) = ^2 Q(l, n, s)N(l, n, s). 



(11) 



4 



As previously shown for square systems, the value of 
X e ff depends on the system's size. We define the ef- 
fective A, X e ff(W, L), as the A for which the analytical 
approximation yields npp(p c ) = 1/2, with p c obtained 
by the numerical simulations. Previous simulations [48[ 
showed that for square systems in the FA model A e // 
does not change much in the region L w 10 2 — 10 5 . Fig- 
ureHJshows that for constant W, the value of X e ff(W, L) 
converges to a finite value A e //(M / , oo), which is rather 
close to X e ff(W, W) at large W. See also Fig. In 
the range of large W and L we see that the value of A 
depends mostly on the model, and not that much on the 
system's size or shape (square or rectangle). Figure [5] 
shows the value of A for very long systems. We see from 
it that for squares it appears that A decreases with the 
width or converges to some value, but for long tunnels 
it is clear that the value of A increases with the width. 
Unless noted otherwise, in the rest of this paper we use 
the value of A taken from the simulations done on large 
squares. 

Figure [6] compares between the fraction of frozen par- 
ticles, Tipp, obtained by the numerical simulations and 
the analytical approximation. From it we see that the ap- 
proximation is roughly acceptable at W — 10, quite good 
even at W = 20, and has an excellent agreement with 
the numerical results at W = 100, even for the hard- wall 
boundary conditions. Also, we note that the approxi- 
mation is better for the periodic boundary conditions, 
and that the width of the transition from an unjammed 
state, where almost all of the particles are unfrozen, to a 
jammed state, where almost all the particles are frozen, 
is narrower with the periodic boundary conditions. 



III. LARGE SYSTEMS, W,L > 1 

In large systems, the transition from jammed (npp ss 
1) to unjammed (npp « 0) occurs in a very narrow region 
of densities. In what follows we find the critical density, 
p Cl at which this transition occurs and the width of the 
transition, Ap. We also show how the critical density 
depends on the shape of the system and not only on its 
size, by considering a system of fixed area and changing 
the aspect ratio W/L. Note that for finite-sized systems 
(and even when L is infinite but W is finite), there is 
no singularity in any physical quantity. Nonetheless we 
use the term critical density since permanently frozen 
particles exist due to the same considerations that govern 
jamming in the thermodynamic limit, where one may 
discuss the notion of critical phenomena [f| 0, H3| • 
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FIG. 4: (color online) The effective A as a function of W/L for 
hard-wall (a) and periodic (b) boundary conditions. In this 
range, A e // does not change much with the length L, and 
converges as L — s* oo to a finite value for each width, W. For 
each system size, A e / / for the KA model is higher than in the 
FA model, and it is higher with hard-wall boundary conditions 
than with periodic boundary conditions, which means that 



n PF > n PF and n 



hard — wall 



> n 



periodic 
PF 



, as expected. At 



L — > oo, A is almost the same for W — 20 and W = 100 but 
different for W = 10. 



use the known critical density for squares of size W x W, 
Eq. ([1]), to find a lower bound on the critical density in 
rectangles of size W x L (W < L) 



v r < 



X 



InW 

From Eq. ([12]) we find that 



(12) 



k = - Win ( 1 - e- 2X/v " ) < We- 2lnW = W' 1 < 1. 

(13) 



A. Critical Density 



Critical Density from Fraction of Frozen Particles 



Since W is very large, this means that the exponent in 
the logarithm is very small, and thus 



k^We~ 2X/v . 



(14) 



We first note that as the system grows larger, the crit- 
ical density grows as well and nears 1. Therefore, we can 



We also note that the critical density is very close to 1 but 



still (pc) W *C 1. Therefore, close to the critical density, 
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FIG. 5: (color online) The effective A as a function of the 
system's width W for hard-wall (a) and periodic boundary 
conditions (b). The plots show the value of A e // at L = 200 VT 
(for W < 100) and L = 5000 W (for W = 200), which is large 
enough to be considered infinite in this range of widths, and 
at L — W. For L S> W the effective A has a minimum 
at W « 50, but there is no drastic change in its value in the 
range 20 < W < 100. Even at W = 200 the relative difference 
from the minimum is 0.08. A e // for long tunnels (L 3> W) 
and squares (L = W) is almost the same at 60 < W < 100. 
We expect that at W — > oo, A will converge to 7r 2 /18 ta 0.54. 



we can use the results from Appendix [X] and approximate 
n PF by 

n PF sa 

1 + cxp [-Lp 2W (kp- 2W + 1)] [kL (kp- 2W + 1) - 1] 



(kp- 2W + l) 2 [1 - cxp (-Lp 2W )] 



(15) 



In very short rectangles, such that Lp 2W <§C 1, we find 
that npF is finite only if kp~ 2W 3> 1 and kL is finite, 
and thus 

n PF ( Lp 2W « 1) « e- kL . (16) 
Solving the equation ripp = 1/2 yields 

2A 



In (WL) ' 



(17) 
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FIG. 6: (color online) Fraction of frozen particles upf vs. 
density p for hard-wall (left) and periodic (right) boundary 
conditions, and for different widths: W = 10 (top), W = 20 
(middle), W = 100 (bottom). Symbols are results of nu- 
merical simulations and continuous line is analytical approx- 
imation with different values of A e //. As the system's size 
increases, the approximation becomes better and the value of 
A e / / becomes dependent mostly on the model, and not on the 
system's size. A was set to X e f f(W, L), as shown in Fig. fj] 



which retrieves the known result, Eq. (JXJ) , for the case 
W = L > 1. 

In very long rectangles, such that Lp 2W ^> 1, we find 
that npp is finite only if kp~ 2W is finite and thus 



n PF (Lp 2W > 1 



1 



(kp~ 2W 

Solving the equation np F = 1/2 yields 



ir 



(18) 



^16XW + ln 2 (W) - In (W) 



AW 



(19) 



This means that the distinction between short and long 
rectangles is whether L {p c ) 2W is greater or lesser than 1. 
Equating Eqs. (IT71) and (IT9l . we find that the crossover 
from short rectangles to long rectangles occurs at 



lnL r = 2Wv r = 



16XW + ln 2 (W) -InW 



(20) 
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2. Critical Density from Typical Section Length 

Another approach for finding the critical density is by 
considering only a typical section of length (/) , where 



vv(WxL) 



(l) = 



Ef =1 Q(0 



(21) 



is the average section length. In this case, similarly to 
what was done on square systems, we need to solve the 
equation 



1 = W (I) e- 2A ^ 



(22) 



Using the expressions in the Appendix, and assuming 
that p w < 1 and L, W > 1, we find that 



-2W 



1 



Lp 



2W 



1 - exp {Lp 2W ) 



(23) 



In the limits of Lp 2W > 1 and Lp 2W < 1, solving Eq. 
(|22|) yields the same results as in Eqs. (fl~7| and (fT9]) . 

This means that when the rectangle's length is shorter 
than (I) it can be considered as consisting of a single 
section, and that for longer rectangles we can consider 
only sections of average length. For this reason short, 
square-like rectangles can be considered to contain only 
one section and the critical density in them depends as 
a first approximation only on the system area WL and 
not on its shape. 

By considering terms of order Lp 2W — L/L c , we find 
that the correction to Eq. (fT7|) is 



2A 



L 



A 



2-ln 2 (2)-ln(4) 



In (WL) L c In 2 (WL) 



ln d (2) 



(24) 



By keeping the ratio L/L c constant, we see from Eq. 
(|2"4")l that the correction becomes less important at larger 
systems. 



3. Alternative Approach on Confinement 

So far we looked at what happens when the width re- 
mains constant and the length increases, and found the 
crossover length L c between the two limiting cases de- 
scribed above. Another way to look at it is by starting 
from a square of size L x L and to generate confinement 
by narrowing its width. Namely, we fix L, and decrease 
W. In this case we find from Eq. (f20| a crossover width 



ln%L) 
4A : 



(25) 



which may be interpreted as the width below which the 
system can be considered infinitely long. To find the 
effect of the confinement on the critical density, we cal- 
culate the ratio between the critical density in squares 




FIG. 7: (color online) Ratio between critical vacancy density 
in rectangles of dimension W x L and critical vacancy density 
in squares of dimension L x L in the KA model with hard- 
wall boundary conditions. Symbols are results of numerical 
simulations and continuous line is the approximations (|17[) 
and (|19|l . The ratio between the critical densities significantly 
differs from 1 only for W < W C (L). W c /L decreases with L, 
and the ratio between the critical densities increases with L. 



and in rectangles. For W > W C (L) or equivalently for 
L < L C (W), this ratio is 



v c (W x L) ^ 
v c (L x L) ~ i + 



1 



21nln(£ c ) ' 
ln£ 



(26) 



2\/AWc 

and for W < W C (L), or L > L C (W), the ratio is 



v c (W x L) 




v c (L x L) 



(27) 



Figure [7] shows the ratio between the critical vacancy 
density in rectangles of size W x L and the critical va- 
cancy density in squares of size L x L. We see from it 
that the asymptotic values in Eqs. (JTTJ) and p^|) agree 
with the numerical results even at W ~ W c . This means 
that systems really may be divided into long and short 
rectangles with a clear distinction between the two types. 

The critical density in Eq. (|19p is a very good ap- 
proximation even for W as small cLS 3 , cLS shown in 
Fig. [H This dependence of the critical density on the 
width of the system in long channels can be measured 
in experiments. The suggested value for the effective 
A is 0.257(FA), 0.275(KA) for hard-wall boundaries and 
0.249(FA), 0.27l(KA) for periodic boundaries, since this 
is its value for systems of infinite length and with a width 
of W « 20 - 100. 



B. Varying Aspect Ratio 

Here we consider a system of fixed volume V = WL, 
and study how changing the aspect ratio between the 
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FIG. 8: (color online) The critical vacancy density, v c , as a 
function of the width for very long systems (L — 200W ^> 
L C (W)) with hard-wall boundaries. The symbols are the re- 
sults from the numerical simulations and the continuous lines 
are the approximation (Eq. (1191) ) with A = A e //(100, 100). 



width and the length, r = W/L, effects the jamming 
transition. The crossover aspect ratio, r c , is the as- 
pect ratio which defines whether the system behaves as a 
square-like system or as a long system. From Eq. (|25[) we 
find that the crossover aspect ratio satisfies the equation 



16Aa 



= In 2 (r c /V) . 



(28) 



If the density is high enough that the system is jammed 
at r = 1, then it is also jammed at any other r < 1. If the 
system is not jammed at r — 1, the density can still be 
the critical density at some aspect ratio smaller than 1. 
This means that as the aspect ratio decreases, the system 
may undergo a jamming transition if the density is below 
the critical density at r = 1 but above the critical density 
at r < 1. 

For example, we now show that our results may be re- 
lated to recent experiments of Bi et al. [I?} • In these ex- 
periments shear stress was applied on a two-dimensional 
system, such that its area and particle density remained 
constant, but the aspect ratio between its length and 
width changed. We consider a system of fixed area 
WL = 10 4 and density p = 0.92 with hard-wall bound- 
aries. The density was chosen such that it is below p c for 
r = 1 and above p c for r <C 1. By changing the aspect 
ratio, the fraction of frozen particles changes from almost 
at r > 0.1 to almost 1 at r < 0.01 as seen in Fig. ® 
The crossover aspect ratio from Eq. (f28|) is r c — 0.093, 
very near the r = 0.1 observed in the numerical results. 
This means that shearing the system causes it to become 
jammed, in agreement with the experimental results pTj . 
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FIG. 9: (color online) Fraction of frozen particles vs aspect ra- 
tio r — W/L at constant area WL = 10 4 and particle density 
p — 0.92 with hard-wall boundaries. Symbols are results of 
numerical simulations for KA (blue squares) and FA (purple 
circles) models. Continuous lines are analytical approxima- 
tions with A e //(100, 100) » 0.271 (KA), 0.252 (FA). 



where 5 is an arbitrary number much smaller than 1. In 
short rectangles [L (p c ) 2W "C 1], setting np^ in Eq. (HI 
equal to S <C 1 and to 1 — S yields 



Ap(short) = p u - pi f 
2A , fhxS- 



ln (WL) 



■In 



2A 



(29) 



where pi^ u are the values of the density at the lower 
(npp = 5) and the upper (npp = 1 — S) bounds. For 
long rectangles [L(p c ) 2W ^> 1], setting npp in Eq. (fT8|) 
equal to 5 and to 1 — S yields 



Ap(long) = p u - pi 



31n<5- 
8W 



3v 2 c hi 5~ 
8T~ 



(30) 



which slightly differs from Eq. (|29l) only in the numerical 
prefactor. 

Figure [10] shows the critical density and width of 
the transition for systems of width W = 100. We 
see, for example for W = 100, that the critical den- 
sity and the transition width indeed converges at L c w 
5000 (iOl), 3000 (FA), in agreement with Eq. O which 
predicts L c = 4606(KA), 3152(FA). For the smaller 
width, W = 20, there is also a convergence, but the fit is 
not as good as for W = 100. The difference between the 
numerical results and the approximations for large W is 
visible because W = 100 is still not large enough for the 
asymptotic limit. 



IV. SYSTEMS WITH SMALL WIDTH 

C. Transition Width 

In systems of small width we can improve our approx- 
We define the width of the transition, Ap as the differ- imation, and even have exact results. In the first two 
ence between the densities for which npp — S and 1—8, subsections below we find an exact result for npp for 
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FIG. 10: (color online) Critical vacancy density (a) and tran- 
sition width (b) with hard-wall boundaries. Symbols are the 
same in both panels. In panel (a), full lines are the approx- 
imations (|17I19I) for large W . In panel (b), full lines are the 
results from the full analytical expression and the dashed lines 
are the large-W approximations. 



L = oo and W = 1 or 2 in both the KA and FA models. 
In the third subsection we improve our approximation for 
systems of width 3 < W < 6 by finding the exact number 
of frozen particles in small sections. 



A. Fraction of Frozen Particles at W = 1 and L = oo 

In systems of width W = 1 we note that the boundary 
conditions are irrelevant. With hard-wall boundaries, a 
particle is blocked from above and below by the bound- 
aries, and with periodic boundaries it is blocked by itself. 



1. KA Model 

For systems with W = 1, we denote by / the number 
of occupied sites in the section and by h the number of 
vacancies in the section. As there are no two adjacent 
occupied sites in the section, in order for a particle to 
be permanently frozen in the KA model, the entire sec- 
tion must be frozen and moreover it must be of the form 



01010.. .01010, i.e. h = f + 1. The number of frozen 
particles in systems with W = 1 is thus 



N«A(p,l)='£v h p h - 1 (h-l) = - 



pv 



h=l 



(i-pvy 



(31) 



and thus, using Eq. (flO)) . the fraction of frozen particles 
is 



n^(p,l) = 2p-p 2 + -^-,. 

(1 - pv) 



(32) 



2. FA Model 

In the FA model, only the particles in the strips are 
frozen for W = 1, and thus 



n^(p,l) = 2p-p' 



(33) 



B. Fraction of Frozen Particles at W — 2 and L = oo 

For systems with W = 2 the boundary conditions are 
important. 

1. KA Model, Hard- Wall Boundaries 

Consider the following two types of patterns: First, a 
pattern with zigzag diagonals of occupied sites, such that 
the other sites are either occupied or not, and second, a 
pattern with a full row of occupied sites, with the sites 
in the other row either occupied or not. In each of these 
two cases, the particles in the main part (the full zigzag 
or full row) are frozen, and in the secondary part they are 
either frozen (if the main part is a zigzag) or unfrozen (if 
the main part is a row). A section with frozen particles 
can be built by dividing it into subsections with their 
main part either a zigzag or a row. Each two of these 
subsections must be divided by a divider, which consists 
of a full column and one particle in each of the adjacent 
columns, one at the top and one at the bottom. The 
following scheme shows this more clearly: 



1010dd000200dd3333 
01010dd2222dd03030 



(34) 



The sites marked with d belong to a divider, and the 
sites marked with a number belong to one of the sub- 
sections. The first subsection is of the zigzag type. The 
second subsection is of the row type. Note that the site 
marked 2 at the top row is not frozen. The third sub- 
section can be of either type. We will consider it to be 
of a zigzag type, since all the particles in the secondary 
part are frozen. In order to simplify the following calcu- 
lations, we will include the rightmost column of a divider 
in the subsection to the right of it. By denoting d as 
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the number of dividers, we note that the number of sub- 
sections is d + 1. We account for the possibility of two 
adjacent dividers by considering subsections of length 0. 
Also, since the left column in the divider is counted in 
it, we need to artificially add the rightmost column in 
the rightmost subsection, since it is not counted in the 
(non-existent) divider to the right of the last subsection. 
Another point to make is that the leftmost column in a 
subsection cannot be full. We denote each subsection by 
the number of vacancies, hi, and the number of occupied 
sites, fi, in the secondary part, such that the main part 



contains hi + fi sites, and by its type, t$ = z,r (zigzag 
or row). A particle in a section built in this way is un- 
frozen only if it is in the secondary part of a subsection 
of row type. We also need to make sure there are no 
two adjacent columns. Also, for each such section, there 
is a mirror configuration with the top and bottom rows 
switched, and so we can count the number of frozen par- 
ticles in one such configuration (say, with the occupied 
site on the leftmost column in the top row) and multiply 
by 2. 

The number of frozen particles in such a section is thus 



d+l 



N 



KAJiw 
PF 



(p,2) = 2 



d=0 



oo hi-S ti ,- 

E E 

i=l ti=z,r ln=0 fi=0 



n e 



v h ip h i+ 2 fi 



fi 



f)V 



d+l 



(35) 



The factor of 2 at the beginning is for the top-bottom symmetry. In the sums, we go over each subsection and check 
how many vacancies and occupied sites there are in the secondary part, where we note that in a subsection of row 
type fi < hi (otherwise we consider it a zigzag type) . The factor of pv before the square brackets is for the rightmost 
column. The sum in the square brackets requires more explanations. First we add the particles in the main part 
(hi + fi)- Next, we say that an occupied site in the secondary part is frozen only if the subsection is of zigzag type. 
The 1 at the end is for the particle in the rightmost column. The final result from evaluating the sums in Eq. ([35]) is 
(see Appendix ICl) 



N 



KA,hw 
PF 



(A 2) = 



2pv (l - p 3 v) [(l - p 3 v) (l + Ap 3 v) + pv (2 - pv + 6p 3 v)] 
[(1 - p 3 v) (1 - 2p 3 v) - pv] 2 



(36) 



and the fraction of frozen particles is 



KA.hw / o\ 



o 3 (2 



P AT KA,hw/ r,\ 
-J N PF (Pi 2 ) 



(37) 



2. FA Model, Hard- Wall Boundaries 



In the FA model with hard-wall boundaries, a section can be at least partially frozen only if all of its subsections 
are of row type, and the frozen particles are only those in the dividers and in the main part. Thus, the number of 
frozen particles is (see Appendix |D|) 



d+l oo hi 



A^ (P ,2)=25>vrn EE'" v 



T FA,hw 

KP,*) = ' 

d=0 i=l lu=Q fi=Q 

The fraction of frozen particles is thus 



pv 



FAM 
l PF 



(p,2)=p 3 (2-p 2 ) + 



3d + X>i+/i) + l 

3=1 



p 8 v (l + 3p 3 v) 
[1 - pv (1 + 2p*)f 



2pv(\ + Zp 3 v) 
[1 - pv (1 + 2p 2 )] 



(38) 



(39) 



3. KA Model, Periodic Boundaries 



In the KA model with periodic boundaries, a section is frozen only if all of its subsections are of zigzag type, and 
thus the number of frozen particles is (see Appendix [0) 



N 



K A,per 
PF 



d-\-l oo hi 



( p,2) = 2£(A) d nEE«v^ 



d=0 



i=l hi=0 fi=0 



pv 



d+l 

3^ + ^(^ + 2/^ + 1 

3 = 1 



2pv (1 + Ap 3 v) 
[l-pw(l + 2p2)] 



r, (40) 
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and the fraction of frozen particles is 

*' per (n.2) = n 3 (2-« 2 l + 

[l-pv(l + 2p 2 ) 



KA,per f , 3/0 2 ^ , p 8 V (l + Ap 3 v) 

pf 0,2) = P (2 - P ) + ; - 9 , l2 - (41) 



«pf 




FIG. 11: (color online) Fraction of frozen particles vs. density 
for the FA (full symbols) and KA (empty symbols) models, 
for hard-wall (squares) and periodic (circles) boundary condi- 
tions, for W — 1,2. Symbols are results of numerical simula- 
tions with L — 200W, continuous line is analytical expression. 
For W = 2 the results for KA with periodic boundaries and 
FA with hard-wall boundaries are almost the same. 

The reason that Eqs. (|3T)|) and (l4"Tj) are very similar is 
that we count almost the same number of particles. The 
only addition to Eq. (|41[) from Eq. is the particles 
in the secondary part. 

4- FA Model, Periodic Boundaries 

In this case, the only frozen particles are in the strips, 
and thus 

4^(p,2) = p 3 (2-p 2 ). (42) 

Figure [TT] shows how perfectly the expressions (l37l) . 
([55]). QSD , and (g2]) fit the numerical results. 

C. Fraction of Frozen Particles at 3 < W < 6 

In principle, we can find exact results also for systems 
with W > 2, but as seen in the previous subsections it 
gets progressively more complicated with increasing W. 
However, we can improve the approximation calculated 
previously. For systems with small width (3 < W < 6) 
we can calculate exactly the average number of frozen 
particles in sections of length I, N(l), by simply count- 
ing all possible configurations. As the number of possi- 
ble configurations rises exponentially with the section's 




FIG. 12: Fraction of frozen particles vs. density for W = 
3, 4, 5, 6 with hard-wall boundaries. Dashed lines are result of 
the approximation of Section [IT] for the KA (higher curve) and 
FA (lower curve) models. Continuous lines are results from 
the approximation of Section IIV CI for the KA (higher curve) 
and FA (lower curve) models. Full (KA) and empty (FA) 
squares are numerical results. Numerical simulations were 
done with L = 200W. Analytical results are for L — oo. As 
W increases, larger sections should be included. For W = 3, 4, 
the current approximation is better than the approximation 
of Section [HI but for W > 5, longer sections are required. 



length and the system's width, we will consider this only 
for sections of size Wl < 40, for which the number of con- 
figurations is 2 40 10 12 , a number which can be handled 
numerically. The frozen particles in the longer sections 
are neglected. 

Figure [T^] shows the fraction of frozen particles ob- 
tained by the numerical results, this improved approxi- 
mation, and the previous approximation. We see that for 
W = 3, it is enough to consider / — 13, but for the wider 
systems we need longer sections. Since the fraction of 
frozen particles in a section decays roughly exponentially 
with the section's length (see Fig. [3]), we need to consider 
only I w 2 (I), where (I) is the average section length, see 
Eq. (j2"3")l . As the average section length depends on the 
density, we can take it at the critical density. We find 
that for such small widths, the average section length is 
6.2 (for W = 3), 6.53 {W = 4), 7.05 (W = 5) and 7.6 (for 
W = 6). This explains why I — 13 is enough for W = 3, 
but Z = 10 is not enough for W — 4. The number of con- 
figurations to scan numerically is 2 2W ^ = 2 39 «5x 10 11 
(W = 3), 2 52 « 4 x 10 15 (W = 4),2 75 « 3x 10 22 (W = 5) 
and 2 96 w 6 x 10 28 (W = 6). 
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V. INTERNAL STRUCTURES 

The main qualitative difference between systems of 
large widths and small widths is the internal structures 
within the sections. In wide systems, the vast majority 
of sections are either almost completely frozen or com- 
pletely unfrozen, while in narrow systems there is a sig- 
nificant number of sections which are partially frozen. 
The reason for this difference is the existence of small 
unfrozen "islands" , which are small regions that do not 
effect their surroundings. For example, a structure of the 



^frozen/unfrozen (^H^ 

1.0 1 — 



form 
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( } q is unfrozen (in both the KA and the FA model) , 

but it does not necessarily cause the entire section to be 
unfrozen. When the system is wide these islands are 
not important, but in a narrow system they are. Figure 
1131 shows the density of completely frozen/unfrozen sec- 
tions in a wide system (W = 40) and in a narrow system 
(W = 7), and Fig. [T4l shows a snapshot of the system, 
highlighting the frozen and unfrozen particles. 

A quantitative way to measure the effect of the bound- 
ary conditions is by noting that with hard- wall boundary 
conditions, the particles in the top and bottom rows have 
a slightly higher probability of being frozen than those in 
the middle rows. This happens because if within a sec- 
tion the top or bottom row is completely full, then all the 
particles in it are frozen, while for the middle rows this 
condition is not sufficient. However, as seen from Fig. [15j 
this difference is very small. Although at first the rela- 
tive difference between the probability of being frozen at 
the edges and at the middle grows with W, this rela- 
tive difference reaches a maximum at W — 12 and then 
decreases for larger W . Hence, we can say that above 
W = 12 the boundary conditions become less important. 



VI. SUMMARY 

In this paper we investigated the effect of confine- 
ment and shape of container on the jamming transi- 
tion. We derived an analytical approximation for the 
fraction of frozen particles in rectangular systems in both 
the Kob- Andersen and Fredrickson- Andersen kinetically- 
constrained models, by dividing the system into indepen- 
dent sections and using the notion of critical droplets, 
which was derived previously for square systems in the 
FA model. The number of these critical droplets is con- 
trolled, in addition to the system's size and particle den- 
sity, by a single parameter A. We showed that the effec- 
tive value of A does not change much when the system's 
length increases, and that its value in rectangular systems 
is approximately the same as in square systems. Also, its 
value in the Kob- Andersen model is higher than in the 
Fredrickson- Andersen model, which means that the frac- 
tion of frozen particles in the KA model is higher than 
in the FA model, as expected by comparing the kinetic 
constraints of the two models. From both the numeri- 
cal simulations and the analytical expressions, we found 
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FIG. 13: (color online) (a) Density of sections which are com- 
pletely frozen or completely unfrozen, hard-wall boundaries. 
The minimum for W = 40 is much lower than for W — 7 
because there is a larger chance of a small, inconsequential 
unfrozen island, (b) The distribution of the fraction of frozen 
particles per section for W = 7 and p = 0.9 and for W = 40 
and p — 0.95 (near the minima in panel (a)). Most sections 
are either completely unfrozen or almost completely frozen. 



that the transition from an unjammed state, where most 
of the particles are free, to a jammed state, where most of 
the particles are frozen, occurs over a very narrow range 
of densities. Using our approximation, we found the criti- 
cal density at which this transition occurs, p c , which con- 
verges to 1 as the system size increases. We also found 
that the width of the transition scales as (1 — p c ) 2 - 

For infinite tunnels, we derived an exact result for very 
narrow systems (widths 1 and 2) for both the KA model 
and the FA model. The technique we used can also be 
applied to wider systems. For infinite systems of width 
3 — 6 we showed that it is enough to explicitly count the 
number of frozen particles in sections shorter than twice 
the typical section length, since the number of frozen 
particles in a section decays exponentially with the sec- 
tion's length. Also, using the analytical approximation 
for general rectangles, we found a simple expression relat- 
ing the critical density with the width of an infinite tun- 
nel, which can be verified in experiments. In particular, 
we found that the critical density in channels decreases 
algebraically with the system's width, 1 — p c ~ 1/ vW, 
much faster than the logarithmic decrease in square sys- 
tems, 1 — p c ~ 1/hxL. These two different expressions 
for the critical density show that the jamming transition 
depends not only on the system size but also on its shape 
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and the relation between the system's width and length. 



(a)W = 100 



(b)W = 7 

FIG. 14: Snapshots of (a) 100 x 2000 system at p = 0.957, (b) 
7 x 140 system at p — 0.85. Both with hard-wall boundaries. 
Black squares are particles which are frozen in both models, 
dark gray are frozen only in the KA model, light gray are 
unfrozen in both models, and white squares are vacancies, 
(a) contains two sections, which are both almost completely 
frozen for KA, but only one of them is frozen for FA. Small 
scattered islands are also visible. In (b) there are 11 sections, 
showing many different behaviors. Zooming in on the pic- 
tures allows to see the details. In print in the 100 x 2000 
system, what appears brightest are the light-gray areas, since 
the individual white sites are too small to be seen. 



Our idea of dividing the system into independent sec- 
tions can also be applied in three-dimensional systems, 
and will be addressed in future work. It will also be inter- 
esting to check the effects of confinement on the behavior 
of other kinetically-constrained models, such as jamming 
percolation models [24], HH , and to employ our approach 
for studying jamming in driven systems, such as granular 
matter flowing in a narrow tube. 
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FIG. 15: (color online) (a) Probability that a frozen particle 
belongs to a certain row at W = 3. At some width-dependent 
density, the probability of a frozen particle to be in the ex- 
treme rows is maximal, (b) The ratio between the maximum 
for the extreme rows and the minimum for the middle row(s) 
as a function of W . This is maximal at W = 12. 



We thank Svilcn Kozhuharov and Vincenzo Vitelli for 
fruitful discussions, and Yariv Kafri and Peter Sollich 
for critical reading of the manuscript. This research was 
supported by the Israel Science Foundation grants No. 
617/12, 1730/12. 



13 



Appendix A: Calculation of the Fraction of Frozen Particles in Eq. (0 



In this section we calculate explicitly the sums in Eq. ([7]). We first note that Q(n, I) is proportional to p n v lw 
such that 



nQ(n, I) = pd p [Q(n, I)] 



(Al) 



where d p [/(/?)] is the partial derivative of /(/?) with respect to p while assuming that v is constant. Only after 
differentiating we use the relation v + p = 1. Hence, npp can be simplified by 



n PF 



En J, m P°P [P^] Q( n > + P mW P d P IQfa 0] e 



-ki 



^[En,I 1 m/* mWr O(".0 

We can now perform the sums over n, such that 

E;, m Pd p [p mW ] Q(0 + P mM V p [Q(0] e- fe ' 



(A2) 



(A3) 



with Q(Z) = Q(l, n). In order to have the derivative in the nominator outside the sum, we artificially change some 
of the p to p' , and after differentiating set p' = p, and thus 



n PF 



pdp 


Ei, m p mW Q(hp') + P ,mW Q(i,p)e- k{ 


pdp 







For ease of calculation, we divide npp into three parts 



npF 



Ai+A 2 
B ' 



(A4) 



(A5) 



where 



Ax = pdp 
A 2 = pdp 
B = pd p 



J2p mW Q(i,p') 

I .III 

Y,p' mW Q(i,p) 

I .III 

"£p mW Q(hp) 



-hi 



(A6) 



and calculate each part separately. 



1. Calculation of Q(l) 
For I — 1 we denote by < / < W the number of particles in the column, such that 

W \ , . -.w w 1 w 

i' • y 

f=o 



Q(l) = £ P f v W ~f { } )=(P + v) W ~P W = l- P w - (A7) 



For the longer sections, we denote by r\ and r 2 the number of particles in the rightmost and leftmost columns, and 
by rii the number of columns with i particles in them. Thus, the density of sections of length I is 

= E E E II E ^+^ 2 — ( £ ) ( I ) p^^ x 



n— ro — nvi/— i— 1 

i 



(A8) 
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where \_x\ is the integer part of x. The first binomial in the second line is the number of ways to arrange nw columns 
among the I — 2 available places such that there are no two adjacent full columns. The upper limit of the sum over 
nw is such because above it the binomial is zero. The second binomial in the second line is the number of ways to 
arrange rij columns among the columns not yet taken by the already-placed columns. Note that we do not sum over 
no, since it must satisfy no = I — 2 — 53i=i n i- 
Summing over ri , f2 yields 



Q(l)= (p + v) w -p w 



n\\r—0 



w 

i 



Wnw -(i-2) (l-nw-1 
n w 



l-2-T w 



(A9) 



We now note that 



W-i^-i:T=l+i n i 

n e 



i=l Ki=0 



A^-V^ 1 n 



j=i+i j | = 
n,: 



w-i 



where go = 1 . The proof for this equation is given in Appendix [B] We can now write 



Q(l)= {p + v) w -p w 



E© 



,W(i-2) 



riw=0 



(p + v) 



W W 



E 



w 

1 + *] I 



Performing the sum over nw yields 



0(0 



(p + «) -p 



tr 



- 1 





X VP 




2 




-?) - 


-1 


-4 











1-p 



w 



y/1 + 2p W - 3p 2W 



w 



where 



/ — nw — 1 

Z — nw — 1 
nw 



w-i 
£<7 

L i=0 



1 + -1- 



(A10) 



-, i-2-r 



py / 



(All) 









2 






h 1 


-4 













I", x w 




2 













(A12) 



i-p w ± \/(i + p w ) -V 



(A13) 



Only in the last step did we use v + p = 1. Using Z = 1 in the general equation for I > 1 yields the same result we 
found in Eq. (|A7I) for Q(l). Note also that Q(l) depends on p and W only via p w . An interesting point to make is 
that 



Q(2)=Q(3)= (1-P W ) 2 . 
Also, we find that Q(7 + 1) > Q(Z) for all values of p , and equality holds only for I = 2. 



(A14) 
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2. The Final Result for upf 



Using Eq. (|A12|) we can now calculate the sums in Eq. (|A4|) . All the sums are such that I + m < L. Also, m is 
greater or equal to 2 except for the case I = L where m = 0. 



a. The Denominator 



The denominator of npp in Eq. (|A4j) is 

B = P d p Y.p mW Q( l >p) 

l.m 

Performing the sum over m yields 

B = pd p 



L-2 L-l 



Q(L,p) + p WL + J2Y,P mW Q( l >P) 



1 = 1 m=2 



Q(L lP ) + p WL + Y i 



L - 2 n 2W n W(L-l + l) 
WL , \ "* P P 



1-P 



W 



Q(1,p) 



Performing the sum over I, differentiating with respect to p and finally setting v = 1 — p yields 
Wp i-2W W p w ^[2-2p^ -p^+L(2-^p w + 3p 2W )} WL 

ri = — — 1 ; — ; ;— Ti h LW p — 



l-p w 

W [x\ + x 1 !] 



{l-pW)(2-3p w y 



2p 2W (l + 3p w )(2-3p w y 
x [2Lp 2W (2 - 3p w ) (1- P W ~ 3p 2W + 3p w+1 ) + (1 + 3p w ) (4p + 2p 3W 

»L1 



(A15) 



(A16) 



jW+l + p 2W+l _|_ p 3W+l 



)] 



+ p w f -Ap 2W W [x^ - x l 
2p 2W (1 _ p W) ( 2 + z p w Z 9p 2W f 

(1 - p w ) (4p + 4p 2W/ + 12p zw - 3&p iW + 54p 5W/ + 12p w+1 - lhp 2W+1 - 54p 3M/+1 + 27p w+1 )] . (A17) 



[2 (L + 1) p 2W/ (2 + 3p w - 9p 2W ) (1 - 3p w + 3p 2W - p w+1 ) 



b. First Part of the Nominator 
The first term in the nominator of Eq. (1A4[) is 



Ai = pd p 



Pdp 



L-2 L-l 



Q(L,p')+p WL + J2T,P mW Q^P) 



1=1 m=2 



= pdp 



L-2 L-l 



y/L 



j2Y,p mW Q^p) 



1=1 m=2 



(A18) 



where in the last transition we note that d p Q(L p') — 0, since here Q is a function of p' and not p. Summing over m 
yields 



Ai = pdp 



L-2 



WL 



£ l 



2W _ p (L-l+l)W 



-Q(1,p') 



1=1 



(A19) 



Summing over I, differentiating with respect to p (not p'), and lastly setting v = 1 — p and p' = p yields 

A W{2-p w ) Wp^ w [L{2-^ + 3p 2W )+3-Ap^} WL 
A l = 1 —/ 1 777-TI 1 77773. h W Lp 



1~P W (1 - p w ) (2 - 3p w f 

W [x\ + x L _] (1 - p w ) (8 - 9p w ) W [x L + - x L _] (1 - p w ) 



7p w - 3p 2W ^ 



2(2- 3p w ) 



2(2- 3p w Y y/i + 2p w - 3p 2W 



(A20) 



1G 



c. Second Part of the Nominator 



The second part of the nominator is 



M = pd p 



E< 



JmW 



Q(l,p)t 



-Ik 



= pdp 



L-2 L-l 



Q(L,p)e- Lk + Y,Y,P 

1=1 rn=2 



ImW 



Q(l,p)e 



-Ik 



= pdp 



Q(L,p)e 



L ~ 2 p/2W _ p/(L-l+l)W 



E 

1=1 



a'W 



Q(l,p)e 



-Ik 



(A21) 



Summing over I, differentiating with respect to p and lastly setting v — 1 — p and p' — p yields 



Nom2 



-kL 



We p 



k n 2W 



£ 2k l p _ p W\ 



1-p 



vl 



We k p W(L+l) (I- pWf + e 2k p W p W^ 



kL p w W [ 



x\ + x^_] 



WL (1 - p) [x£ + x_] 

2(l + 3p w ) T 2(l-p^) 2 (l + 3/) H ') A 

{e 3fe (1 - p w ) [(1 + 3p w ) (1 - 4p w + 4p 2W - p w+1 ) — L (l — p w — 4p 2W 

[yl-yt]+e 2k (l-p w ) 2 [2 + 4p w -6p™ 



!l2 



6p3w + p w + l _ 3p . 



2W+1 



)] 



LI 



p w_ 3p 2W + pW+i)] [ p w yl _ yl] 



-e ik (1 - p w ) [2 + ip w - 6p 2W -L(2 + p w - 6p 2W + 3p w+1 )] [yj - p w y 2 ] + 
+e 5k [1 + 2p w - 3p 2W - L (1 - 3p 2W + 2p w+1 )} [y\ - p 2W y 2 ] } - 
e -kL W [ X L _ X L] ]^ p w (l-p)- L (l + p- 2p w ) (1 + 3p w )} 



2(1 + 3p w ) y/1 + 2p w - 3p 
Wp w [ic£ - x*l] 



,2W 



2 (1 - p w ) (1 + 3p w ) VI + 2p w - 3p 



„2W 



{_ e 2fc(l_ p ^) [2 + 6p w_ 6p 



2W 



6p 



3k 



1 - P W ) [1 



10// 



6p 6 



+ 4p 



2W+1 



L (1 + 3p w ) (l~p w ~ P 2W + P w+1 )] [p w yl - y 2 ] 



3p 



2W+1 



-L(l + 3p w ) (l-2p w + p w ^)] [yl-y 2 ] 



+e 4k [2 + 6p w - 10p 2W + 6p 3W - Ap w+1 - L (l + 3p w ) (2 - 3p w + p w+1 )] [y{ - p w y 2 ] - 



-e 5k [l + Ap w -3p 



2W 



2p 



w+i 



L (l + 2p w — 3p 2W )] [yl- P 2W yl]}, 



(A22) 



where 



1 



Vi 

V2 



(e k + p w ) (1 - p w ) - e 2k ' 
1 

(e k + 1) (1 - p w ) - e 2k p w ' 



(A23) 



Appendix B: Proof for Eq. (TaTo)) 



Here we prove that 



AT 

W-l N -l^j = i + l ™3 

n e w 

i=l m=o 



Hi 



W-l 

1 + E 9i 



N 



(Bl) 



We do this by induction on W. For W = 2 it holds because 

N-0 



ni=0 



E [9irl N -° ) = [i+, m] N 



r>i 



(B2) 
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For W + 1 we first sum over n\ 



i=l n i= o 



[i+.^n e 



i=2 ra»=0 

We now define /i,- such that 















.1 + 91. 





i=2 n i= 



9i+i 



1 + .91 ' 



and rewrite the sum as 



w- 1 "-z^j^i+i "j 



n E &r ( N ~ <* ) = [i + 9l f n E ( " " E r i+1 

i—1 n ,=0 ^ 1 ' i=l n;=0 



Since the sum is now only up to - 1 we know what it is 

n e w 

i=l n 4 =0 

as required. 



'-ij=i 

m 





W-l ' 


N 






)=[l+9if 


i=l 




1 - 


i=i 



(B3) 



(B4) 



(B5) 



(B6) 



Appendix C: Derivation of Eq. (|36p 

The sum in Eq. (I35[) over each subsection is identical and independent of the others, thus we can transform Eq. 
3 to 



A^(p,2) = 2p*£(pV) 



h-8 tl : 



d+1 



(oo h-8 t ,r 
E E E <V +2/ ( ; 
t=z,rfc=0 /=0 

Calculating the sums over /, /i and t yields 



t= z ,/ /i=0 /=o 
d 



h-5, r 



h EE E '-''^ -''( ; ) / + 

, t=z,r h=0 /=0 ^ ' 



(CI) 



N^ hw (p, 2) = 2pv £ (p 3 ^ [(3d + 1) C 2 d+1 (p) + (d + 1) C 2 d (p) d (p)] , 
pu [2 + 3p 2 - 6p 5 u - 2p 4 u 2 - 2pV + 3pV] 



C 2 (P) 



(1 - p 3 v) 2 [1 - pv (1 + p 2 
1 + pu (l - p 2 ) 



(l-p 3 i;) [1 - py (1 + p 2 )] 



Summing over d yields Eq. (|36l) 



(C2) 
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Appendix D: Derivation of Eq. (|38[) 

We rewrite some of the p in Eq. (|38)l as p' , so that Npp can be simplified to 

-ifh 



oo d~\- 1 oo /i ■ 

N*$**{p,2) = 2J2 II E E vWpr<pf"+'«+*+ 

d=0 i=l fci=0 / 4 =0 

oo d+l OO /li ; 



d+1 
3=1 



d=0 i=l fe 4 =0 /i=0 



As before, since each of the d + l sums over hi and /j are independent, we can write Npp as 

d+l 

OO /l / 7 \ 

7V£^(p, 2) = 2p'cV ^ « d + V 3d+1 



d=0 



/!=0 /=0 



Calculating the sum over h and / yields 

oo 

Np ? P hw (p,2) = 2p'd p ,J2v d+1 P /3d+1 



1 



1 - vp' (1 + pp') 



d+l 



Calculating the sum over d, differentiating with respect to p' , and finally setting p' = p, yields 

2pv (l + 3p 3 v) 



N™< hw (p,2) 



[1 - pv (1 + 2p2)f 



(Dl) 



(D2) 



(D3) 



(D4) 



Appendix E: Derivation of Eq. (1401) 



We rewrite Eq. (|4U1) as 



d+l OO /li 



iV 



#^( P ,2) = 2P 9 P ^ (Ar n E E «V i+2/i ( /; 



d=Q i=l h i= fi=0 

This is exactly Eq. (|D1|) with p' = p, and therefore we can use Eq. (|E)3|) 



7V^ CT (p,2) = 2 P 9 p g^+V 3d+1 



d=0 



1 



n d+l 



l-vp{l + p 2 ) 



Summing over d and differentiating with respect to p yields 



[l"P«(l + 2p 2 )] 2 



(El) 



(E2) 
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